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Applying arguments based on the operator product expansion for a three-point correla- 
tor and relying on quark-hadron duality, we derive an expression for the skewed (non- 
forward) parton distribution in the pion in the case of a zero-skewedness parameter, 
y^ll{X;t). We expect that our result is relevant for moderately large momentum transfers 
1 < t < 10 GeV^ In addition, we construct a purely phenomenological factorized model 
for the same quantity in close analogy to Radyushkin's model, originally proposed for 
skewed distributions of quarks in the nucleon. Though the quark-hadron duality approach 
supports theoretically the factorized model, the two models exhibit a different behavior 
in the parton momentum fraction X at any fixed t. The relevant process to distinguish 
between the two options seems to be the WAGS off the pion that measures (to leading t/s- 
order) the inverse moment < > of the skewed distribution. Even after the inclusion 
of the first order kinematic t/s-corrections, the predictions for the cross section ^(s,t) at 
c.m.s. scattering angles = 30° and 90° differ by factors 3.5-3.9 and 2.9-7.5, respectively, 
so that a discrimination appears possible. 



I. INTRODUCTORY REMARKS 

Recently, Radyushkin |^ has argued that at moderately large momentum transfer t — {1 — 10) GeV'^, 
hadronic form factors and wide-angle Compton scattering (WAGS) amplitudes are dominated by a soft 
mechanism corresponding to an overlap of soft wave functions. This analysis was performed in terms 
of universal nonforward parton densities J-'{X;t), which accumulate the soft contribution in the WAGS 
case. These densities are obtained in the C = limit of the nonforward parton distributions (NFPDs), 
J^^{X;t), introduced in and represent the simplest hybrid distribution interpolating between usual 

parton densities f^{x) and hadronic form factors. A simple factorized model for T{X;t) in the proton 
was constructed in ||l[ using the well-known Gliick-Reya-Vogt parameterization |0] for | proton ( 3^) ^-nd 
assuming a Gaussian dependence on the transverse momentum k±_ of the effective two-body soft light-cone 
wave function ^(a;, k±) of the proton. 



*0n leave of absence from the Institute for Nuclear Research and Nuclear Energy, 1784 Sofia, Bulgaria. 
^The NFPDs are similar to but not coinciding with the off- forward parton distributions (OFPDs) introduced by 
Ji in §11; cf. the discussion in 



In this paper, we perform a similar analysis relating to the case of the nonforward parton densities in 
the pion. We show that the factorized model can be approximately justified within the Operator Product 
Expansion (OPE) in conjunction with QCD sum rules [||. 

As a first step towards a complete QCD sum-rule analysis, we explore in this paper the so-called local 
quark-hadron duality approximation (simply abbreviated as LD in the following) ||^-[r^ that was suc- 
cessfully applied to estimate various nonperturbative characteristics, like hadron masses, leptonic widths, 
electromagnetic form factors of hadrons, etc. pO|-p^. Note that our approach has the advantage of being 
simultaneously gauge and Lorentz invariant from the outset. 

Within this framework, we are able to obtain a compact expression for the NFPDs, like !FQ=o{X;t), in 
the form of an overlap integral of the Drell-Yan- West-type [^,|9| , which involves an effective two-body 
soft pion wave function ^'(a;, kj^), introduced earlier in We compare the LD expression for Ti^=o{X; t) 
with the corresponding factorized ansatz and comment on the reliability of the LD approximation. 

Finally, an estimate of the leading contribution to WAGS off the pion is made for moderately large 
scattering angles (in the cm. frame). We also discuss possible corrections of order of 0{t/s) at leading 
twist and show that new nonperturbative quantities should be introduced in addition to the skewed parton 
densities. These new quantities are particular y-moments of the same underlying double distributions 
F{x,y;t). 

The analysis of these moments, as well as an extension of our approach to the nucleon case, will be done 
elsewhere. 

The paper is organized as follows. In Sect. II we establish our definitions of double distributions (DDs) 
for the pion case. We derive their symmetry properties and relations to the corresponding NFPDs. In 
Sect. Ill the operator product expansion approach for the skewed distributions in the pion is discussed in 
connection with local duality. The prediction for the distribution J-'i^^q{X; t) following from this approach, 
is given in Sect. IV. The consistency of the LD result with the general sum rule for the pion form factor 
is demonstrated. Sect. V deals with the modeling of the skewed distributions, using a phenomenological 
parameterization for the valence quark distribution in the pion. We demonstrate that the WAGS process 
off the pion is very sensitive to the model for Tc^^Q{X]t), especially to its small X behavior. In Sect. VI 
we discuss finite 0{t/s) corrections to the handbag contribution of the skewed distributions. Finally, in 
Sect. VII we further discuss our results and draw our conclusions. 

II. DOUBLE DISTRIBUTIONS IN THE PION; DEFINITIONS AND SUM RULES 

First, we define double distributions in the pion in close analogy to the nucleon case [|^,^ in terms of 
a nonforward matrix element of a bilocal quark- antiquark operator on the light cone using, however, the 
most general decomposition (see, also [pO|) 

{tt{p') I i^iOh^EiO, z; A)^{z) \ ^(p)) U^^o = 

{P + P')^j^ I e(x + y<l)(e-"(^'^)-^^(")F^I^(x,y;i)-e"(^'^)-^5(-)^^l^(3,^y.i)) dxdy 

+ {p- p')^ J^' e{x + y< 1) (e-"(f^)-*«(-)G^|,(x, y; t) - e^^^^^^-^y^"-^^ G^^,{x, y; t)) dx dy (1) 

+ ^^^^f^f^ e{x + y<l) (e-"(^'^)-^^(")Z^|^(x,2/;t)-e"(P^)-^^('-^)Z^|Jx,y;0) dxdy, 

where r = p — p' is the momentum transfer (r^ = — t < 0) and E{Q,z]A) = P exp(ig dXfj_A^{x)) 

is the path-ordered gauge string factor in the fundamental representation {A^ = X]a=i ^a^/j)-EI ^^-(0) 
corresponds to operator structures of leading and next-to-leading twists (twist 2 and 3, respectively). 

Matrix elements of this type appear in the perturbative QGD analysis of deeply virtual Gompton scatter- 
ing (DVCS) processes QJ^. One can actually prove in all orders of perturbation theory that the large-Q^ 



'Note that to leading- twist accuracy, the definition of the matrix element in Eq.(|l|) is independent of the choice 
of the contour connecting the quark fields tlj{0),ip{z) (cf., e.g., |21[). 
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asymptotics of the scattering amplitude can be represented in a factorizcd form with the short-distance 
part calculated perturbatively (see, e.g., |p2|,^). The dynamics of large distances, which is mainly nonper- 
turbative, is in turn accumulated in the matrix elements of the type given by Eq.(^. 

The z^-term in Eq.(Q) is of higher twist-3, that is, it will produce a power-suppressed contribution 
(~ l/Q^) in the DVCS amplitude. Another two-body twist-3 part is produced by the axial- vector matrix 
element 

{irip') I V:;(0)7p 75i?(0, z; A)^(z) | ^(p))U.^o = 



ifj. 

r-l 



"'0 



+ y<l)(e-"(^'^)-2^(")A^|,(a:,y;i)-e"(^'^)-*^('-^)A^|.(a;,y;t)) dx dy . (2) 



Because we neglect in this paper power-suppressed corrections, our main focus will be on the F-, G- 
distributions. 

The parameterization of the nonforward matrix element in terms of DDs is natural and can be es- 
tablished, at least in perturbation theory, at any order of Us Just as in the case of deep-inelastic 
scattering, DDs in Eq.(|l|) have a parton interpretation: y] t) is the amplitude to find an active quark 
in the pion with momentum fractions x and y of the initial (hadron) momentum p and the momentum 
transfer r, respectively. 

Due to the conservation of the local vector current, the DDs G^\^{x, y; t) obey a sum rule. Indeed, for 
z = only the first term in Eq.(|l|) should survive (r^ ^ 0) and this implies 

0{x + y<l)(^G^^^{x,y;t)-G^^^{x,y;t)')dxdy^Q. (3) 
Moreover, there are stronger sum rules: 



0{x + y<l)G^^^^^{x,y;t)dxdy^O. (4) 

Indeed, using TP-invariance and complex conjugation of the matrix elements in Eqs.(|l|),(|^), one can easily 
show that the DDs introduced above are real- valued functions and should obey the symmetry relations 

F, p\^{x,l-x-y;t) ^ F^\^{x,y;t) 

G, p\^ix,l-x-y;t) = -G^\„{x,y;t) (5) 
Z^\^{x,l-x-y;t) ^ Z^i^{x,y;t) 

A^l^{x, I - X - y;t) = A^\^{x, y; t) , 

which are a generalization of the "Miinich symmetry" relations ||2^ in the case of the pion. Hence, the 
sum rules encoded in Eq.(^ are an obvious consequence of that symmetry for the DD G^\t^(x, y; t). 

Within the generalized Bjorken limit for DVCS, we have the relations Q^,pq' ^ t, rrv^ and r|| = (p [| 
where C = Q^/2pq' coincides with the Bjorken variable. Considering C as an external parameter, one can 
introduce the NFPDs ||^,^ T(^{X;t), with X = x + (^y being the total fraction of the momentum of the 
active parton, to read 

J^p'^iX; t) - e{X > / F^\^{X ~ 2/C, y; t) dy + e{X < C) / F^\^{X - yC, y; t) dy . (6) 
Jo Jo 

Analogously, one can also define the "forward invisible" NFPDs gf'^iX^t), zf^iX-.t), Af^'iX-t). 
Then the leading contribution of the handbag diagrams in Fig. |l] can be obtained in the form 



"(C - 2) (t^^{x- t) + Tf{x- 1)) - c {gfix- 1) + gf{X; t)] 



(7) 



nl 

dX 



X-C + ie X-ie 
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FIG. 1. Diagrams contributing to the DVCS amplitude. The blobs at the bottom corre- 
spond to DDs. 



Another process, described with the same handbag diagrams, is the wide-angle Compton scattering 
(WACS), recently reexamined in papers ||l|j2^ for the proton case. Now the initial photon is also real 
(Q^ = 0), but t = is large enough to ensure the light-cone dominance (see [Q for a discussion of other 
contributions with subleading 0{t/s) behavior). The contribution of large distances in this case will be 
described by the same skewed distributions at C = 0, notably 

Tfil{X;t) = T^\^{X;t). (8) 

Of course, in this case one should use Eq.(|^) with some care because finite i-corrections to the hard part 
may become important [^,0. These t-corrections, which are within the leading-twist approximation, 
are analogous to the target-mass corrections in DIS which have led to the Nachtmann-Georgi-Politzer 
^-scaling ]2|f|l. 

As it was shown in refs. I2l0], DDs play a key role in describing those processes in which nonforward ma- 



trix elements are involved. In fact, many properties of the skewed distributions, like polynomiality 28 1, 
symmetry properties [ p3| , p8| , etc., can be simply established using the integral representations of the type 
ofEq.®. 

However, DDs seem to have a more complicated structure. In fact, the question of possible singularities 
of the DDs is still open. A pure perturbative analysis seems to yield DDs without singularities However, 
DDs are by definition nonperturbative objects and one may expect such singularities to appear in the real 
world. In the work of |28[ , t-channel meson-exchange contributions were considered and found to produce 
(5- function type singularities. In the case of pion DDs a dynamical mechanism was found [po[ , based on 
the effective chiral model which follows from the instanton vacuum of QCD | |29| , |30t . 

Independently of whether or not these singularities really appear in the DDs, it is clear that the skewed 
distributions, which are certain integrals of the DDs, are more smoothly-behaved functions. Thus, they 
seem to be more appropriate for modeling. 



III. OPE APPROACH TO THE SKEWED DISTRIBUTIONS 



In this and the subsequent sections, we are going to study the skewed distribution in the pion, J-^^^ {X; t), 
within an approach which is based on QCD sum rules with particular emphasis being placed on the 
case C = . As in the proton case [^jsj , one can derive a sum rule connecting the charge pion form factor 
with a certain integral of the NFPD (the zeroth moment), namely. 



1 ^ 



Tf^{X- 



t)-Tp''{X;t) 



dX 



(9) 



where 60 is the electric charge of the active "?/;" -quark (see also Q). Indeed, Eq.(^ follows immediately 
from Eq.(|l|) after taking the limit z = and incorporating the sum rule for the DD, termed G in Eq.(|^). 
It should be emphasized, however, that for fixed ^ 0, the physical domains for the t-variable on the 
left and the right hand side of Eq.(||) are different. In fact, for the DVCS process, the limit t — > is 
unreachable fB2l, i.e., 
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(10) 



On the other hand, in the forward hmit {( = 0, t = 0), a reductfon formula holds, (see and we 

have 

T^il{X;t = 0) = f.^\^{X). (11) 

Both, the form factor, as well as the parton distribution in the middle region of X, were thoroughly 
investigated within the QCD sum- rule approach. In fact, for the pion form factor F^{t) it was shown that 
in the region of momentum transfers, t > 1 GeV^, the so-called Feynman mechanism is capable to 
reproduce the experimental data |jl^,^ without recourse to the hard part. 

In this paper we shall adopt a similar philosophy and use to derive the skewed distribution in the pion the 
concept of local quark-hadron duality ||]- [|l^ that was successfully applied to the calculation of various 
nonperturbative characteristics, like hadron masses, leptonic widths, electromagnetic form factors, etc. 
(for some recent applications, we refer to ||l5| , [l^ , |l6| , |l7| ) . 

As usual in the QCD sum rule approach, let us consider the three-point amplitude 



{p,p'-z) = i^ jd^x j d^y e-^P^e^P'y {^()\t [£ {x)mi^m.z-.A)u{z) jl{y)]\Q) , (12) 



Tin 

where j^{x) — d(x)j5jau{x) is the axial current with a nonzero projection on the pion state, so that 

(0 |ja(0)| 7T+{p)) = tUpa. , U - 133 MeV, (13) 

and z denotes a light-like coordinate (z^ = 0). 

The correlator defined by Eq.(|T^) will be considered in the Euclidean region for p,p',r = p — p' . It 
gives contributions to different invariant form factors with tensor structures proportional to Pa{p + p')fj.p'[3, 
Pa'Tfj.p'p, gafsTfj,, etc. For the C = case we shall project on a light-like direction via nanpn^, where v? = 0, 
p+ = (np), and r+ = {nr) = 0. The advantage of this projector is that it projects out the leading structure 
in the infinite momentum frame (IMF), where — > oo, with r±_ fixed. This structure is also most directly 
related to the one analyzed before in by employing the QCD sum-rule method to calculate the pion 

form factor. 

Picking out the invariant amplitude of the leading structure, we have 

where the ellipsis denotes polynomials in p^,p'^. The perturbative contribution to i?ti(p^,p'^, t; z) (which 
is the leading term of an OPE expansion in the deeply Euclidean region of the momentum invariants) can 

be written in the same form as (^ with the obvious change pu^-^^ pS*^"^*- Due to asymptotic freedom, 

for large s and s', p^^^ (s, s') ~ /9S^'^*(s, s'). However, for small s,s', the two spectral densities differ 

drastically from each other. Indeed, pu^^^ contains the pion double (5-function term 

pZ{s, s', t, z) = 2^V^ {s ~ ml) 5^+^ {s' - ml) 

(e-'^(P-).F^l;(X; t) - e'^(P^).Fg,(X; t)) dX , (15) 

whereas, in contrast, f^^^{s, s') is a smooth function for any finite order of perturbation theory. 

Under the proviso of the local duality (LD) assumption, one has that pj^(s, s') is dual to p^^^{s, s') in 
an appropriate duality interval, so that 



- / pr\s,s')dsds' ^- / pl[s,s')dsds'. 
Jo Jo Jo Jo 



(16) 



Here the duality interval sq corresponds to the effective threshold of the higher excited states and the 
"continuum" in the channels with the quantum numbers of the axial current. 
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The LD relation ( |16[ ) is very natural within the QCD sum-rule approach g]. In fact, the effective 
threshold Sq is fixed by the ratio of the nonperturbativc power corrections (the condensate contributions) 
relative to the (leading) perturbative term in the OPE for the correlator (cf. Eq.(12)). In what follows, 
we shall use the value 



So — '■'0 

which follows from the LD prescription for the correlator of two axial currents: 



^Mp)^i e'^P^ (o|T{jf (x)4(0)}|o) 



(17) 



(18) 



For the experimental value of we have Sg^ = 0.67 GeV^. This value is very close to the standard 
one, 0.7 GeV'^, that has been extracted from the direct QCD sum-rule approach for the 2-point 

correlator, Eq.(^8|), in the classical work of Ref. in determining the pion decay constant /tt. 

The same duality interval was also obtained in the QCD sum-rule analysis of Refs. |^,^ of the charged 
pion form factor at moderate momentum transfers t « 1 — 3 GeV^. However, it was observed in that 
for higher values of i, the relative contribution of the condensate (power) corrections increase, ensuing 
an increase of the extracted parameter so as well. This situation corresponds to the so-called infrared 



regime [^,^, i.e., to the kinematical regime, where one of the quarks carries most of the momentum of 
the initial hadron (current). In this regime the underlying OPE series becomes badly convergent and should 
be resumed in some way, e.g., by introducing nonlocal condensates {q{0)q{x)), {G{0)G{x)), {q{0)G{x)q{y)) , 
etc. [^^. Adopting a reasonable model for the nonlocal condensates, it was shown that the form 
factor i^7r(i), extracted from such an improved approach, can describe the data up to i ~ 10 GeV^ in 
compliance with previous rough estimations ]39|JTl[ |, according to which the asymptotically leading hard- 
scattering contribution starts to become important beyond t ^ 10 GeV^. As discussed more fully in ]35|] , 

the corresponding threshold s^^ ^'^^ < * < 10 GeV ) g^^j^ ^ type of QCD sum-rule analysis was found to 
have approximately the standard value, quoted above. 



IV. LOCAL DUALITY PREDICTIONS FOR THE SKEWED DISTRIBUTIONS IN THE PION 

The one-loop contribution to the double spectral density can most easily be calculated using light-cone 
variables in a frame where the initial momentum p has no transverse components, i.e., 

P= - ^, 0^1 (19) 

and the momentum transfer has no "plus" component (cf. |[T^): 

P' 

Here, s, s' are the invariant masses in the channels with pion quantum numbers, X is the total fraction of 
the longitudinal momentum carried by the quark entering the composite vertex, and kj^ is its transverse 
momentum. 

Applying the Cutco sky rules (cf. Fig. ||a), one can obtain for the double spectral density of the "i?u"- 
correlator (recall Eq.(|l^)) in leading order of as 



0, r--4-^^±ii, f^l , fc={Ap+, fc-, k^} . (20) 



I 



\ 



k± - ?±X 

Yx 



2 ^ 



d^kx . (21) 



For the contribution at hand, the spectral constraint < X < 1 reflects the positivity of the energy of 
the struck and spectator quark, respectively. In fact, this constraint has a more general nature, discussed, 
for instance, in P,p[. 
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FIG. 2. Typical Cutcosky cuts (dashed lines) for perturbative diagrams in the OPE for 
the 3-point correlator of two hadron currents, involving a composite operator (crossed oval 
at the top). The left and right graphs correspond to the 0(1) and 0(as) contributions, 
respectively. 



Substituting Eq.(pl|) into the LD relation, provided by Eq.(|T^) one can extract (by taking the Fourier 
transform) the corresponding skewed distribution 



/ 



e so 



XX 



e 



So 



v^X 



\ 



XX 



(22) 



Analogously, one can deduce that, in leading order, the following relations hold 

■^c5),fi|7r(-'^; ^) = •^c5),d|7r(-'^; ^) = 0- 



(23) 



Eq.(|2^) can be recast in the form of an overlap integral of effective two-body soft light-cone wave 
functions of the incoming and outgoing pion: 



■L4„|^(X;i)= / ^ 



/.LD 



167r3 



(24) 



The explicit form of the effective wave function ^jj^^jx, kx) can be obtained, in particular, from the 



evaluation of the 2-point correlator, Eq. (|T^) , to read 



V'e^D(X,k^) = 



V6 



e Lso - 



XX 



(25) 



Thus, the LD prescription reproduces an expression of the Drell-Yan-West type In a recent 

paper , the classical Drell-Yan-West formula was used as a starting point for developing a phenomeno- 
logical model for the form factors and skewed distribution of the proton. 

Carrying out the integration over the transverse momentum k^ , we obtain from Eq. (M) 



e X 



t 



4so + t 



soXX 



/ Xt 



/ Xt 



Xt 



(26) 



Note that the 8- function in Eq. ( p6| ) is due to the corresponding abrupt behavior of the effective pion wave 
function ( p5| ) dictated by local duality at leading order. It ensures that the corresponding X^^-moment 
oi T^^{X]t), which enters the leading handbag expression for the WAGS, is indeed finite (cf. Eq.(0)). 
We expect that at next-to-leading (NLO) order in as (cf. Fig. ||b) of the LD-approach this property 
will remain intact. In fact, such a property is required by the general considerations of the factorization 
theorem M^. 



7 



For small momentum transfers t, the form of the quark-hadron duality will change (see the discussion 
bellow). However, in a more realistic model for the effective wave function and/or skewed distribution, 
this property should be preserved as well. 

Taking the zero-order X-moment of ^ul^{X;t) and incorporating the relations (p3|), one can reproduce 
to leading-order accuracy the general sum rule Eq.(p|) (at C = 0) 



Jo 



(27) 



where 



^ ^'^"'"(l + 4Vt)^/- ^''^ 

Actually, the explicit LD prediction for the pion form factor was obtained earlier in |P,[T^. 

Thus, the concept of local duality provides us with nontrivial dynamical information about the form 
of the pion wave function/skewed distribution, as well as about the behavior of the pion form factor at 
moderately large momentum transfers t > 0.6 GeV^. As we shall see in the next section, this form- factor 
prediction seems to be supported by the existing experimental data in this momentum region. Moreover, 
due to the Ward identity, connecting the 3-point function Eq. ([l^) and the 2-point function, Eq.(|l8|), the 
property F^(i = 0) = 1 is also preserved]^ 

On the other hand, it is well known that in the region of small momentum transfer t, the quark-hadron 
duality is more complicated pO|,p|,p 15 1. Thus, one should not overestimate the accuracy of Eq.(p7|) in the 



region of i < sq. In fact, the derivative of Eq.(27) at zero momentum transfer is infinite. The reason is that 
in the kinematical region t ^ |-P^ | i \p'^ \ ; 0^*^ has to include additional terms in the OPE corresponding 
to the situation in which the current with small momentum transfer is placed at a large distance. This 
leads to the notion of bilocal power corrections |^,^. Taking them into account extends the validity of 
the theoretical QCD sum-rule prediction for the pion electromagnetic form factor to the whole momentum 
transfer region i = (0 — 3) GeV^, providing, in particular, the correct value of the pion charge radius 

It is important to note that the same reasoning can be applied to a more complicated object, namely, 
the skewed distribution Tu\Tr{X; t). In fact, the t = limit of Eq.([2^) gives 

.FLD(x. i ^ 0) - /LD(X) = QXX. (29) 

The parton distribution f^^{X) is normalized to unity and hence respects the same Ward identity, men- 
tioned above. Moreover, it turns out that its form coincides with the asymptotic (leading- twist) distribution 
amplitude of the pion ip'^{x). Thus, as a consequence, the (naive) quark-hadron duality procedure fails 
to reproduce a reasonable valence parton density in the pion. Actually, to reproduce the small i-behavior 
of the skewed distribution T^^{X;t) within a QCD sum rule approach is a rather complicated problem. 
Both, the notions of bilocal power corrections of leading twist , as well as non-local condensates [btIJs^ 
should be introduced. We shall address this interesting problem in a forthcoming pubhcation. 

However, for large momentum transfers t >\ GeV^, in analogy to the form factor calculation, one may 
expect that the LD result, given by Eq.(pq), should work. 



V. MODELING THE SKEWED DISTRIBUTIONS 

In order to get some experience of how reliable the LD strategy is, and to estimate how large the 
deviations in the region of small momentum transfers t may be, we shall derive in this section a factorized 
model for the skewed parton distribution in the pion at ^ = 0. 



^This property was shown to be fulfilled in a complete QCD sum-rule analysis for the pion form factor, performed 
in H. 
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We shall argue that the factorized ansatz for ; t) , following the approach of Ref. [Q, is of the 

form 



which makes it apparent that it automatically satisfies the general "reduction relation" (Ill| 
that is the valence u-quark distribution in the pion 



(30) 



Note 



The specific functional dependence of the exponential in (30) on t- and X 
within the OPE approach. In fact, if one replaces into the overlap formula 
function, i/;LD(X,ki), given by Eq. (p5|) , by the popular Gaussian ansatz, proposed in 



Gaussian 



(X,ki) = <i>(X)e-'^i/2A'^* 



can be formally justified 
2^ ) the abrupt LD wave 

1,11, 

(31) 



one also arrives at such an exponential dependence. Another hint at an exponential dependence can be 
traced back to the double Borel transform of the 3-point correlator (^2|) when employing the OPE. In 
fact, the t-dependence of the perturbative term (first diagram on the rhs of Fig. |^) and that of the term 
involving a (vector) nonlocal quark condensate (inserted into the bottom line of the second diagram on 
the rhs of the same figure) is described by one and the same function, namely. 





f + 




, + 



FIG. 3. Operator product expansion of the 3-point correlator of two pion currents involving a composite 
operator. 



$(t,X) = e-*^/(^^?+*-^')^ . (32) 

Here Ml and M| are Borel parameters, corresponding to and p'^ , respectively, and X is the momentum 
fraction, flowing through the upper lines. Other terms of the nonlocal OPE - not displayed in the figure 
- are more complicated, but numerically their f-dependence is similar. 

A Gaussian distribution form for the transverse momentum should not be surprising. It is in line with the 
Borel transformation technique, as it was recently shown in p^ , where effective Gaussian wave functions 
were derived. 

The dimensional parameter A in (^o|) should be determined in correspondence with some averaged values 
of the Borel parameters {Mf 2) for which the underlying sum rule is saturated. Thus, one may expect that 

A^ « (Aff 2) « 2 • 0.7 GeV^ = 1.4 GeV^ , (33) 

where 0.7 GeV^ is the characteristic scale for the two-point pion correlator and its 3-point counterpart is 
twice larger |35| . 

On the other hand, A can be determined on a purely phenomenological ground by adopting the Gliick- 
Reya-Schienbein (GRSch) parameterization for the valence quark distribution in the pion | p7| : 

(a;) = 0.5645 a;°-5°4(l + 0.153 VS)(l-a;)°-^*^ , (34) 

which refers to a rather low normalization point ^^o ~ '^•^^ GeV^. According to our considerations, we 
shall use in our modeling procedure a parton distribution evolved to a more appropriate scale: /i^ ~ 1 GeV^. 
Taking into account the leading order LD estimates for ^^lo(A;i) (recall Eq.(|2^)), one is tempted to 
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FIG. 4. Predictions for tF^(t) based on the LD model, Eq.(Eq), (solid line) 
and the factorized GRSch-ansatz^ Eqs.(p^-(p5[) (dashed linejT^ The experi- 
mental data are taken from 14511 



assume that the valence d distribution is f^^^{x) ~ fuiiri^) and that the sea quark contributions are 
practically negligible. (This assumption seems indeed to be supported by experiment | ^ , ^ ). 

To fix the parameter A, we use sum-rule relation (|^), related to the electromagnetic pion form factor, 
and a properly weighted sum of quark and antiquark DDs. Taking into account (E3f) and e„ — = 1, one 
can write 



Factorized 



7 .T--ii[7r;Factorized / ^ ,\ 
dA y-^^p (A ; t) 



( 



tX \ 



(35) 



which makes the dependence on the parameter A explicit. The best agreement between our factorized 
(and OPE-inspired) model (pq) and the experimental data WSMGl in the region of intermediate momentum 



transfer 1 GeV^ < ^ < 10 GeV^, is realized for A§ w 1.7 GeV^ (see Fig. y) which is not far away from the 
QCD sum-rule inspired value A^ — 1.4 GeV^, given above. 

Now, we are able to compare the two models in more detail. The results are presented in Fig. ^ for 
three different values of the momentum transfer: < = 1, 3, 10 GeV^. One observes from this figure a quite 



t = 1 GeV 




t = 3 GeV^ 




0.2 0.4 0.6 O.S 

X 



t = 10 GeV^ 




0.2 0.4 0.6 0. 

X 



FIG. 5. Predictions for the skewed distributions in the pion with C = 0, obtained with local duality (solid lines) 
and the factorized ansatz in conjunction with the GRSch parameterization (dotted lines), and the (naive) LD 
parton distributions, Eq.(b!l), (dashed lines). 



different behavior over the momentum fraction X (at fixed momentum transfer t) for the two displayed 
models, though the general tendency is that for higher t the distributions tend to shift their weight towards 
the upper end-point of the interval < X < 1. For the sake of completeness we also present in the same 
figure, a factorized model based on the naive LD distribution function, Eq.(p9|). As one may expect, this 
option gives a somewhat "smoother" behavior relative to the LD curve. 

By construction, the zero-order X-moment following from the factorized ansatz for 3-'^}^Q{X;t), (cf. 
Eq.(35)), fits the data for the pion form factor rather well. On the other hand, also the LD prediction for 
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the pion form factor (Eq.(p8|)) complies with the data quite well. Thus, in order to distinguish between 
the two models, one should look for other physical observables, which are more sensitive to the form of 
the distribution. 

In principle, the form of the true distribution at given t can be reconstructed having recourse to higher 
moments: < X'^ > {N > 0). In addition, the small X-behavior of J^f}^Q{X;t) is sensitive to the inverse 
moments < X"^ >. Remarkably, a possibility to measure the < X^^ > is offered by the WACS process 
in the pion case. Indeed, as we have mentioned above, the leading part of the handbag contribution to the 
WACS amplitude is proportional to the integral 

< X-' >EE ^ , (36) 

as it can be seen from Eq.(0) at C = 0- 

In the case of unpolarized initial photons, the differential cross section reads 

^i.s;t) = ^{el + eiy{Rl,{t)+Rl,it)f , (37) 

where we are still neglecting finite i-corrections. As we see from this equation, the inverse moment X~^, 
or equivalently i?"i(i), enters squared in the cross section formula. Thus, one may expect an improved 
sensitivity of the WACS process to different models of the skewed distribution in the pion (as well as also 
in other hadrons). 

In Fig. II we show the differential cross sections of WACS off a pion as functions of the squared center-of- 
mass energy s for different values of the scattering angle i? which is hold fixed. One infers from this figure 
that the two models presented above yield cross sections which differ from each other by approximately 
one order of magnitude. 

We remark that estimates for the pion Compton scattering were also presented on phenomenological 
grounds in pOj . 




FIG. 6. Leading order predictions for the differential cross section of WACS off 
a pion, plotted against the squared center-of-mass energy s for different values of 
the (fixed) scattering angle ■& (in the center-of-mass reference system) . The solid 
line shows the result derived from local duality; the dashed line that following 
from the GRSch-based model. 

As we have mentioned in the Introduction, in order to ensure light-cone dominance, one should consider 
a Compton process at momentum transfers t^l GcV^. Thus, the t/s ratio may not be small. We recall 
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in this context that t ~ sin^(z9/2)s in the center-of-mass reference system, we are using, and d is the c.m.s. 
scattering angle. 

In such a situation, to make a more reliable comparison with the experimental data, one should actually 
take into account finite t/s corrections. This task will be undertaken in the next section, where we consider 
possible corrections of order of 0{t/s) at the leading- twist level. 

VI. FINITE 0{t/s) CORRECTIONS WITHIN THE HANDBAG CONTRIBUTION 

In order to derive the 0{t/s) corrections at the leading- twist level, one can start from the contribution 
of the handbag diagram in the coordinate representation (see, ||l5[| ). Taking into account the parameter- 
izations of the leading twist non-forward matrix elements through DDs (F'^ (x,y;t), {x,y;t), etc.), via 
Eq.([^), and performing a Fourier transformation, one arrives at the following expression for the Compton 
scattering amplitude 



T^.{p,q',r)^i J2 4 J dxj dyeil>x + y)S^p,p[ip + p')p(^F^ix,y;t) + F^ix,y;t)^ + 

+ [p- p')p (C'^ix, y- 1) + G"^"(x, y- i)) } {Sp{qs) - Sp{qu)) + "A-term" + "Z-term" , (38) 



where S^p^p = Qtipgi^p - g^ugpp + Qfj^pgi^p and Sfj{qs), Spiqu) are, respectively, the hard quark propagators 
for the s and u channel handbag diagrams: 

Spiq)^^, qs = q + xp + yr, qu = q - xp + yr . (39) 

q 

For real Compton scattering, i.e., q^ = 0, the denominators of the quark propagators can be written in 
the form 

ql — xs — y{l — x — y) — xx and q^ = xu — r^y{l — x — y) ~ xx (40) 

and respect the "Munich symmetry" , Eq. (|^) . 

In principle, it would be legitimate to retain in the calculation (see als o [|l5|l) 0{m'^)-, 0(r^)-terms as 
well, in analogy to target- mass effects in DIS that have lead to ^-scaling [ p5|j2(^ ] . Here, in the pion case, 
the rn^-terms are not significant and can be neglected compared to s, t, and u. 

Suppose for a moment that we neglect the i-corrections to the hard quark propagator. Then, we 
immediately reproduce the leading handbag contribution Tp^{p^q' ,r), Eq.(^, with the leading tensor 
structure —gfn, -\- -^{Pfiq'i^ + Pvq' p) depending on the two Sudakov 4- vectors p, (?'. The appearance 

of the non-forward (skewed) distributions JF^'''(X;i), CJ^'^(X;t), defined through the DDs F^|^(x, y; t), 
G^l^(x, y; t), in Eq.(|^), is in accordance with the fact that in the formal t limit, the denominators of 
the quark propagators depend on the combination X = x + C^y only, the latter being the total momentum 
fraction of the active parton. For C, — Q, we have 

T^\^X;t)= dyF^\JX,y;t) , etc. (41) 
Jq 

In general, when t 0, the scattering amplitude Tp^{p,q' ^r) should depend on three independent 
4- vectors which may be chosen to be p,q',r±, or p,q',r, etc., and will also include non-leading tensor 
structures.^ The denominators will also acquire an additional nontrivial dependence on the x- and, 
y-fr actions. 

It is in practice more convenient to consider directly the (i/s)-corrections to the cross section of WAGS. 
For the case of unpolarized initial photons, and summing over the polarizations of the final one, we obtain^] 



**In fact |rx| ~ Vi- 

The SU{2) symmetry relations (cf. Eq.(p3)) have been taken into account. 
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da , 



t r 

s 



+ Ki,it)y {Ri,{t) + Ri,{t)) (Rut) + Rm) 

i 



where R\{t) is the corresponding first X-moment 

< >= R'iit) = [ dXXT\ 
Jo 



(42) 



(43) 



The first term in Eq.(^) reproduces the leading result for the cross section, i.e., Eq.(p7|). It is worth 
remarking that the leading t/s-corrections are expressed only through the moments of the skewed parton 
distribution T^]^^{X]t). In contrast, the NLO t^/s^-corrections include all distributions introduced in 
Section II. 

Let us now introduce the y-moments of the DDs 



nX 

T]^^{X-t) = dyy^F^\,{X,y-t) 



(44) 



and analogously for the other DDs, namely, G^|^, ^^i^ri ^i/>|7r- Then the skewed distribution J^^^'^{X;t) 
becomes simply the zeroth-order y-moment of the corresponding DD, i^^|^(X, y; t). Moreover, one can 
show that the t^/s^-corrections can be expressed by means of the X-moments of these new distributions. 
Obviously, the higher y- moments, Eq.(^), cannot be expressed in terms of the simplest skewed distribution 

JP'^l'^(X;t). Thus, one is forced to introduce new skewed distributions, ^Ij^y^^^ {X] t), though these distribu- 



tions refer to the same DD. The generalization of Eq.(| 
but we shall skip it here for the sake of brevity. 



to non-zero skewedness, C 7^ 0, is straightforward 
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FIG. 7. Local duality (LD) predictions for the differential cross section of 
WAGS off a pion with the first order kinematic t/s-corrections included as a 
function of the squared center-of-mass energy s for selected values of the fixed 
scattering angle (in center-of-mass reference system) i9. 
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The evaluation of the different y-monicnts is a separate task which we shall not address in this paper. 
Thus, we restrict our analysis here to the first order i/s-correction to the cross section, Eq.(|2|). The 



evaluation of the moments Ri{t) is straightforward by virtue of Eqs.(p6D, (30). 

Fig. shows the plot of our LD predictions for the differential cross section of WACS off a pion as a 
function of the squared center-of-mass energy s for different values of the fixed scattering angle (in the 
center-of-mass reference system) d. We expect that the prediction is reliable for i ^ 1 GeV^, as it was 
explained above. Notice here that the condition t^l GeV^ for the kinematics with fixed scattering angle 

transforms into the condition s > 1/ sin^(i?/2) GeV'^. For the minimal (30°) and maximal (90°) values 
of the scattering angle d this means s ^ 15 GeV^ and s > 2 GeV^, respectively. 

As it turns out, the relative magnitude of the i/s-correction to the LD-predicted skewed distribution 
J^^°(X;i) is « 4% for 1} = 30° and rises to 10-28% for i9 90°. On the other hand, in the case of the 
factorized ansatz, (X ; t) , the corrections are stronger (5% and 25-42%, respectively). 

Still, even after taking into account the t/s-corrections, the predictions for the cross section in the LD 
case, evaluated a.t ■& — 30° (90°) is 3.5-3.9 (2.9-7.5) times smaller than with the corresponding factorized 



ansatz. Hence, forthcoming experiments at TJNAF |49| may be able to discriminate between these two 



models. The most dramatic difference of the two models appears, however, in the kinematical region 



d = 90° and for s > 10 GeV^ 



VII. DISCUSSION AND CONCLUSIONS 

Let us start our discussion with some comments concerning the energy-momentum sum rule for the pion 
parton distributions, which wc have rewritten in the form: 

dXX (j^"l'^(X;t) + J^'^>(X;i)+ J^"l^(X;t)+ J^''l^(X;i)+ J^sk(X;t)) It^o = 1 • (45) 

Here we have neglected the sea quark contribution of other flavors and J^3^^{X;t) is the C = version of 
the gluon skewed distribution in the pion that can be defined in an analogous way to the quark case (cf., 
e.g., §). 

With the leading order LD result (meaning zeroth order in a^), i.e., Eqs.(|26|), (|23|), taken in the t = 
limit, one observes that the sum rule, Eq.(^) is already saturated with the quark contribution alone; the 
sea and gluon contributions being exactly zero in this approximation. 

In this way we demonstrate the consistency of the LD approach with the energy-momentum sum rule, 
Eq.(45). It is interesting to mention that the energy-momentum sum rule is saturated by the quark 
contribution alone also within the effective chiral model of Ref . |2^ . Actually, their approach is based on 
a low-energy effective action derived from the instanton vacuum (see, p9|j30|] and references therein) and 
states that the gluon distribution should be parametrically small, ~ (p/i?), where p « ^ MeV~^ is the 
average instanton size and R is the average distance between instantonsj^ 

We would like to emphasize once more that the LD approach described in Section III is, strictly speak- 
ing, not applicable to the region of small momentum transfers t < 0.6 GeV^ and should therefore be 
modified. For instance, using our factorized model for J^'^^'' {X]t), Eq.Q, which is based on the GRSch- 
parameterization of the quark distributions in the pion, one finds (at pr :± 1 GeV'^) for the quark contri- 
bution a value of ^ 0.54 to the energy-momentum sum rule with some room for the gluons left over. Note 
also that the contribution of the sea amounts only to ~ 7% of the total result. 

In the region of applicability of our LD formula Eq.(|2^), t ^ 1 GeV^, we expect that the skewed 
distributions of the sea and the gluons should be suppressed. Indeed, within the LD approach they will 
both appear first at the level of a^-corrections. 

Another point worth to be discussed concerns the estimation of the intrinsic transverse momentum of 
the pion. With the explicit form of the effective two-body wave function ^p^'^ {X,'k±), given by Eq.(|2£ 
one immediately obtains 



'■"'•The quark skewed distribution in the pion has been studied in the instanton vacuum, but in a somewhat different 
approach, also in Ref. |feo|. 
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LD 



LD 

-2^ w (260 MeVf . (46) 



On the other hand, employing the factorized ansatz, Eq.(pO|), one gets for the valence quark distribution 

,^ „, Factorized „ /"^ _ „ 

(k^2)^ = / XXf^\^{X) dX « (320 MeYf , (47) 

Jo 

which is in moderate agreement with the LD estimate. In both cases a reasonable value of the intrinsic 
average transverse momentum of the pion is obtained. This value turns out to be somewhat smaller than 
that proposed by KroU and collaborators |Q. The reason for this discrepancy may be traced back to 
the fact that the LD two-body wave function, Eq. (p5|) , is an effective one and therefore includes not only 
the lowest particle-number (quark-antiquark) Fock state, but also an infinite tower of Fock states with 
additional soft gluons (cf. the discussion given in [^Hj and also [^). 

As a result, our ijj^^ {X,'kj^) can account for the full normalization condition 



dX(P]^i_ 
167r3 



V-. (x,ki)|' = 1 , (48) 



whereas the two-body wave function of Ref. |Q contributes approximately 1/4 of it. 

In conclusion, we have provided evidence that the local duality approach for the skewed distribution 
JP^iyX;i) in the pion seems to support both, theoretically and phcnomenologically, the so-called factorized 
model for the same quantity (proposed in ||l| for the case of the proton). However, the two presented models 
show distinct behaviors in the longitudinal momentum fraction X for any fixed momentum transfer t. 

Furthermore, we have shown in this paper that measuring the Wide Angle Compton Scattering off the 
pion will provide us with a very sensitive tool to test the form of the skewed distribution and, in particular, 
to discriminate between the two models, proposed in this paper, and against others. We emphasize that 
our analysis was performed by taking into account the finite t/s-corrections to the cross section of WAGS. 
Thus, we expect that our results may be relevant for an experimental check in the energy region of 
s ~ 2 15 GeV^ and a scattering angle d ~ 30° 90° in the c.m.s. which seems accessible to the TJNAF 



machine |49 



In a forthcoming publication we plan to generalize the LD approach to the case of a non-zero skewedness 
parameter C,. 
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